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Although cosmic microwave background (CMB) anisotropy data alone cannot constrain simul- 
taneously the spatial curvature and the equation of state of dark energy, CMB data provide a 
valuable addition to other experimental results. However computing a full CMB power spectrum 
with a Boltzmann code is quite slow; for instance if we want to work with many dark energy and/or 
modified gravity models, or would like to optimize experiments where many different configurations 
need to be tested, it is possible to adopt a quicker and more efficient approach. 

In this paper we consider the compression of the projected Planck CMB data into four parameters, 
R (scaled distance to last scattering surface), l a (angular scale of sound horizon at last scattering), 
Qbh 2 (baryon density fraction) and n s (powerlaw index of primordial matter power spectrum), all 
of which can be computed quickly. We show that, although this compression loses information 
compared to the full likelihood, such information loss becomes negligible when more data is added. 
We also demonstrate that the method can be used for scalar field dark energy independently of 
the parametrisation of the equation of state, and discuss how this method should be used for other 
kinds of dark energy models. 

PACS numbers: 98.80.Es,98.80.-k,98.80. Jk 



I. INTRODUCTION 

Dark energy model buildingcontinues to be an active 
area of research 0, 0, |, i, MM S i, i, E3- While cur- 
rent data remain consistent with a cosmological constant 
explanation for dark energy, other possibilities are not 
yet ruled out, especially if theoretical motivation can be 
found to tighten their predictions about the data [ill ]. 
New theoretical ideas thus may bolster support in favour 
of an exotic component of matter or a modification of 
gravity beyond some length scale. 

On the observational front, recognizing the need for 
better data, many future dark energy surveys have been 
proposed, classified by the Dark Energy Task Force as 
stage III and stage IV experiments [l2|. The realisable 
constraints from these surveys depend sensitively on the 
external or prior information that will be available in the 
future. A crucial external data set will come from the 
Planck satellite, which will place strong constraints on a 
range of cosmological parameters. It is therefore impor- 
tant to include this data for forecasts and optimisations 
of instrument performance for the stage III and IV dark 
energy surveys. This in turn requires a rapid way to eval- 
uate the predicted Planck likelihood, preferably without 
the necessity to run a Boltzmann code. 

Some of us have shown that the information from 
the WMAP CMB experiment [H[ can be effectively and 
simply incorporated into a likelihood analysis of Type 
la supernovae (SN la) and baryon acoustic oscillation 
(BAO) data by including in the likelihood a term in- 
volving WMAP constraints on the CMB shift parame- 



ter (R), the angular scale of the sound horizon at last 
scattering (Z a ), and the baryon density flbh 2 , in Gaus- 
sian form together with their full covariance matrix [lij . 
The idea being that the calculation of full CMB spec- 
tra at each parameter point can be avoided, so that a 
Markov Chain Monte Carlo (MCMC) analysis proceeds 
very quickly. The merit lies in the method being inde- 
pendent of the dark energy model used as long as only 
background (or homogeneous) quantities are varied. 

In this paper we extend the method to projected 
Planck data, which is significantly more accurate than 
WMAP data. We derive and test this simple prescrip- 
tion, compare it to a full likelihood analysis of simulated 
Planck data, and conclude that when such a Planck prior 
is combined with future dark energy surveys useful com- 
plementary information from the CMB is retained and 
there is no significant information loss. Hence the pre- 
scription remains an effective way to incorporate con- 
straints from Planck (or Planck priors) in the analyses of 
data from future dark energy surveys. 



II. COMPONENTS OF THE PROPOSED 
PLANCK LIKELIHOOD 

Let us first introduce the parameters that we are 
proposing to use as an effective summary of the infor- 
mation contained in a CMB spectrum: 



I? 



y/fim-ffo r(zcMB), la = TTv(zc M b) /r s (zc Mb), 



(1) 
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where r(z) is the comoving distance from the observer 
to redshift z, and t s {zcmb) is the comoving size of the 
sound-horizon at decoupling. We give the details of the 
formulae used in appendix lAl 

In this scheme, l a describes the peak location through 
the angular diameter distance to decoupling and the size 
of the sound horizon at that time. If the geometry 
changes, either due to non-zero curvature or due to a 
different equation of state of dark energy, l a changes in 
the same way as the peak structure. R encodes similar 
information, but in addition contains the matter den- 
sity which is connected with the peak height. In a given 
class of models (for example, quintessence dark energy), 
these parameters are "observables" relating to the shape 
of the observed CMB spectrum, and constraints on them 
remain the same independent of (the prescription for) 
the equation of state of the dark energy. Furthermore, 
R and l a are very well constrained by WMAP and even 
better by Planck and their likelihoods are almost per- 
fectly Gaussian (remaining so under different treatments 
of dark energy), so that a Gaussian likelihood term to- 
gether with the corresponding covariance matrix retains 
almost all of the information on these derived parame- 
ters. With curvature held fixed, an even simpler set up 
using just R sufficed and has been used by many authors, 
including [lH. fl^|. 

As a caveat we note that if some assumptions regarding 
the evolution of perturbations are changed, then the cor- 
responding R and l a constraints and covariance matrix 
will need to be recalculated under each such hypothesis, 
for instance if massive neutrinos were to be included, or 
even if tensors were included in the analysis [28]. Fur- 
ther R as defined in Eq. fTJ) can be badly constrained and 
quite useless if the dark energy clusters as well, e.g. if it 
has a low sound speed, as in the model discussed in [3(| . 
However, as discussed further below we checked that our 
constraints are valid at least for scalar-field dark energy 
models, independent of the parametrisation of w(z). 

In addition to R and l a we use the baryon density flbh 2 , 
and optionally the spectral index of the scalar perturba- 
tions n s , as these are strongly correlated with R and l a , 
which means that we will lose information if we do not 
include these correlations. 



III. SIMULATED DATA 

Our simulation and treatment of Planck data is as in 
[r^ ]. We include the temperature and polarization (TT, 
TE, and EE) spectra from three temperature channels 
with specification similar to the HFI channels of fre- 
quency 100 GHz, 143 GHz, and 217 GHz, and one 143 
GHz polarization channel, following the current Planck 
documentation, 1 . The full likelihood is constructed as- 
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suming a sky coverage of 0.8. We choose a fiducial model 
close to the WMAP best fit LCDM model: n b h 2 = 0.022, 
n m h 2 = 0.127, h = 0.73, Q k = 0, w = -1, and w a = 0. 

For the Baryon Acoustic Oscillation part, we use the 
experimental configuration outlined in the DETF report 
[12| . A Stage III spectroscopic experiment would cover 
2000 square degrees with a redshift range of 0.5 < 1.3, 
divided into 4 equally sized redshift bins, plus 300 square 
degrees with 2.3 < z < 3.3. The experiment would ob- 
tain the spectra of 10 7 galaxies. This survey will mea- 
sure the oscillations in the galaxy power spectrum, in 
the tangential direction (measuring r(z)), and the ra- 
dial direction (measuring dr(z)/dz = c/H(z), providing 
a direct measurement of the Hubble parameter) . To esti- 
mate the accuracy with which the radial and tangential 
oscillations can be measured, we apply these survey pa- 
rameters to the fitting formulae described in [l8[ . These 
fitting formula only consider the accuracy with which the 
oscillations themselves can be measured, returning no 
information about the accuracy of the power spectrum 
measurement (as is done in e.g. [lj||2(J. This is because 
the number of possible parameters contributing to the 
nature of the matter power spectrum, such as running 
of the spectral index, massive neutrinos, and non-linear 
bias, make this calculation very assumption dependent. 
In contrast, the positions of the oscillations is very robust 
with regard to these extra considerations. 

For the Supernovae, we use a Stage III spectroscopic 
survey as described in the DETF report [1J. We as- 
sume a scaled-up version of the SNLS survey with 2000 
supernovae in the range 0.1 < z < 1, with a further 500 
supernovae at low redshift. The dispersion in observed 
magnitude is the sum in quadrature of a fixed <jd = 0.12 
with a second piece <r m , which is fixed at 0.02 up to 
z = 0.4 but then increases up to 0.03 by z — 1. 



IV. ANALYSIS 

The full set of constraints on all parameters includ- 
ing R and l a are determined through an MCMC based 
likelihood analysis [2l| of simulated Planck data. Planck 
will provide much tighter constraints on parameters, and 
its posterior will be significantly better localized in pa- 
rameter space than that of WMAP. The shape of the 
posterior (i.e. parameter correlations) is also found to 
be quite different from WMAP's (further justifying the 
exercise of determining the best way to incorporate con- 
straints from Planck separately from WMAP). While R 
and l a were almost uncorrelated for WMAP data, this 
is no longer the case for Planck. Tables 1 and 2 show 
the estimated values and the covariance matrix for i?, 
l ai flbh 2 and n s . We have included n s here because it 
is found to have a correlation with R and l a and a dif- 
ferent consideration of BAO data in the future, utilising 
the full shape of the matter power spectrum, might re- 
quire the inclusion of n s as a parameter. In the analysis 
presented in this paper, given the conservative treatment 
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FIG. 1: This figure shows the projected 68% and 95% Planck 
constraints on R and l a obtained assuming that dark energy 
were due to a cosmological constant (flat ACDM, dashed con- 
tours), a wo,w a model (shaded contours) and a kink model 
(solid contours), as described in the text. 



of the BAO signal, the inclusion of n s does not have a 
noticeable impact. 

The first point to consider and re-test with Planck 
data is whether the constraints on R, l a , fl^h 2 and 
n Sl and their corresponding covariancc matrix are inde- 
pendent of the dark energy prescription used. In [3| 
we tested this for WMAP data for a cosmological con- 
stant, constant w and wo-w a models of dark energy, 
with and without curvature. Here we test it again for a 
flat model with a cosmological constant, and the wo,w a 
model and the kink model for dark energy, both with 
curvature, and for Planck quality data. In the Wo,w a 
model wx{z) = wq + w a (l — a) [24| which corresponds 
to X{z) = a -3(i+™o+™„) e 3™a(a-i)_ In t he kink model 
the equation of state parameter wx is described by its 
value today, wo, its asymptotic value at high redshift, 
w m , as well by two more parameters giving the location 
and speed of the transition from w m to wo [25j ] . In this 
case the energy density is derived through a numerical 
integration of the continuity equation. We found that 
there is no significant difference in the constraints on R, 
l a , Slfc/i 2 and n s obtained using these different models. 
See Figure 1. 

Let us now test for the amount of information on pa- 
rameters relevant to dark energy that is lost by consider- 
ing a likelihood based on R, l a , flbh 2 and n s rather than 
the full CMB spectra. Figure 2 shows WQ,w a contours ob- 
tained from a full likelihood analysis of Planck simulated 
data (shaded contours) against contours reconstructed 
from the R, l a , Vl^ti 2 and n s likelihood (solid curves). 
We find that even in this limited 2D view there is sig- 
nificant information loss: The shaded contours from the 
full likelihood cover significantly less area than the open 
contours from the simpler likelihood. Due to the strong 
degeneracies which leave Planck basically unable to con- 
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FIG. 2: This figure shows the Planck projected constraints on 
wo,w a obtained using a full likelihood analysis of simulated 
Planck data (shaded contours) and a simpler and quicker like- 
lihood analysis based on R, l a , 0.bh 2 and n s (solid contours). 
Information is thus lost by the simplified analysis. 



strain cosmological parameters relating to dark energy 
and curvature on its own, the resulting contours depend 
strongly on the priors used. 

It may be useful to note that given the R, l a , flbh 2 and 
n s likelihood, one can implement a full likelihood anal- 
ysis under different dark energy models more efficiently 
using Hamiltonian Monte Carlo. In this method the R, 
l ai flf,h 2 and n s likelihood is used as a guide to or an 
approximation of the true likelihood surface, but at each 
accepted point the likelihood is weighted using the full 
CMB spectra. We discuss this procedure in more detail 
in appendix [Cl 

However the question remains, whether there is still 
an information loss when using the R, l a , flbh 2 and 
n s likelihood from Planck when analysing SN la and/or 
BAO data, as compared to the full Boltzmann analysis 
of Planck data, which is much more time consuming and 
so limits our ability to consider many varied dark energy 
models. To address this we compared the outputs from 
two analyses. Firstly we performed a full MCMC run, 
i.e. including the full Planck likelihood and likelihoods 
from simulated stage III SN la and BAO surveys. Sec- 
ondly we perfomed a MCMC analysis using the R, l a , 
Q,bh 2 and n s likelihood from Planck together with the 
SN la and BAO likelihood. Fig. [3] shows the constraints 
obtained in each case. We conclude that there is effec- 
tively no information loss in using the R, l a , Vlbh 2 and n s 
likelihood, in conjuction with the likelihood from a better 
SN la and/or BAO experiment, and this condensed data 
analysis proceeds much faster than an analysis involving 
the full CMB likelihood. 

Another way to include a Planck prior in forecasting 
constraints from a future dark energy experiment is to do 
it via a Planck Fisher matrix. Consider the above like- 
lihood analysis in comparison to a Fisher matrix treat- 
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FIG. 3: This figure shows 2D confidence contours obtained 
using a full likelihood analysis of simulated Planck data in 
conjunction with stage III SN la and BAO data (shaded con- 
tours), contours obtained using the simplified R, l a , fib/i 2 
and n„ based likelihood analysis of Planck data together with 
stage III SN la and BAO data (solid contours) and finally 
a Fisher matrix treatment of all data (dashed contours) as 
described further in the text. 



ment. The Planck Fisher matrix was obtained from the 
Planck covariance matrix of (i?, l a , flbh 2 , ns), with the 
appropriate parameter transformations for compatibility 
with the SN la and BAO Fisher matrices. See Appendix 
B for a description of how the Planck Fisher matrix was 
obtained and Table IV for the resulting Planck Fisher 
matrix. Constraints on dark energy parameters obtained 
in this way are also shown in FigO Because of the nearly 
unconstrained directions, the pure Planck Fisher matrix 
cannot be inverted, as the range of eigenvalues is larger 
than its precision. This can be rectified with weak priors 
on the parameters (in which case diagonal entries in the 
inverse of the Fisher matrix will reflect those priors), or 
by adding more data. Figure 3 shows that the Fisher 
matrix is valid in spite of its formal problems: the error 
contours for Planck + SN-Ia + BAO data agree very well 
with the others. 



V. CONCLUSIONS 

We have found that a Gaussian likelihood based on R, 
l ai fl(,h 2 and n s effectively summarizes the information 
in Planck that is relevant for an analyses of data from 
SN la and BAO experiments for dark energy parameters 
under different dark energy models. Therefore a Planck 
prior can be included in this manner. When used in 
conjunction with other data that are more sensitive to 
dark energy such a treatment of Planck data results in 
no information loss as compared to a full analysis, while 
being much faster. 

We provide the full i?, l a , Qbh 2 , n s covariance matrix 
that is required to define such a likelihood from Planck. 
We also provide a Planck Fisher matrix for people who 
prefer to use the Fisher matrix route to forecasting con- 
straints for a future experiment. Using such a Planck 
prior we have obtained the type of constraints that may 
be expected from a stage III SN la and BAO survey. Of 
course the prescription can also be used once data from 
all these experiments have actually been obtained (ie. the 
prescription is not just for forecasting). 

In the above analysis we found that it was not strictly 
necessary to include n s given our conservative treatment 
of BAO data. A fuller treatment of BAO data such as one 
that included the shape of the matter power spectrum 
rather than the transverse and line of sight distances to 
the redshifts of the BAO survey deduced from the BAO 
scales in the corresponding directions, would require the 
primordial power spectrum parameters including n s to 
be considered a variable in the BAO part of the analysis. 
For this reason we have included n s in our prescription, 
and marginalized over it in our results. 

While this work was in progress 28] considered a like- 
lihood analysis involving the locations of the peaks and 
troughs in the CMB spectrum observed by WMAP to 
constrain dark energy parameters in combination with 
recent BAO data. This offers another way to include 
information from the CMB in a likelihood analysis of 
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BAO and SN la data. It involves fitting formulae for 
the locations of the extrema presented in [29(. Fitting 
formulae have been derived to account for certain pre- 
recombination effects that via the early ISW effect can 
effect the position of the first peak relative to the higher 
peaks. In our formalism we would have to recompute 
the -R, l a , Cibh 2 and n s constraints for each new pre- 
recombination scenario, such as involving a non-zero neu- 
trino mass, involving tensors and/or the running of the 
scalar spectral index, or else include these parameters 
in the covariance matrix. On the other hand, our ap- 
proach is arguably simpler to implement, and at least as 
accurate within its domain of applicability (since it ad- 
ditionally uses R as an effective measure of peak height). 
In their approach too new fitting formulae would have to 
be derived new effects in different scenarios that haven't 
been considered in the past. 



Acknowledgements 



is given by 



r s {z G MB) 



c s dt 
a 

O-CMB 



= cH n 



cm 



dz — — 

ZCMB E ( Z ) 

da 



3(1 + i? 6 a) a A E 2 {z) 



where a is the cosmic scale factor, clcaib = 1/(1 + 
zcmb), and a i E 2 {z) = £l rad + £l m a + £l k a 2 + Q x X(z)a i . 
The radiation density is computed using the Stefan- 
Boltzmann formula from the CMB temperature, assum- 
ing 3.04 families of massless neutrini. The sound speed 

is c s — 3(1 + Rb a), with R^a = 3/?b/(4p 7 ), Rb = 
31500fVi 2 (T CMS /2.7K)- 4 . 2 
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APPENDIX A: DETAILED DESCRIPTION AND 
FORMULAE 

The Planck satellite will deliver data of such a high 
quality that even small changes in parameters like the 
CMB temperature can have a significant impact. For this 
reason we summarise here the relevant formulae used in 
this paper. Generally they are those used by CAMB. 

The comoving distance to a redshift z is given by 

r{z) = cH^ 1 lOfcl-^shmpfcl 1 /^)] (Al) 
r(zH JaW) 1 E(z)=H(z)/H 

where il k = —k/H^ with k denoting the curvature con- 
stant, and sinn(a;) = sin(x), x, sinh(a;) for Q k < 0, 
fife = 0, and f2fc > respectively, and 

e(z) = [a„(i + ^) 3 + a-ad(i + 2) 4 (A2) 
+n k {i + z) 2 + n x x( z )] 1/2 (A3) 

with Qx = 1 — — f^rad — and the dark energy 
density function X(z) = px(z)/ px(Q)- 

We calculate the distance to decouplingm, zcmb, via 
the fitting formula in [H. CAMB [23[ uses the same 
fitting formula. We note that simply using a constant for 
zcmb results in a shift in the inferred values of the CMB 
shift parameters at levels of precision corresponding to 
Planck. The comoving sound horizon at recombination 



APPENDIX B: FISHER MATRIX APPROACH 



The Fisher matrix, F a p , for a set of parameters p can 
be derived from the Fisher matrix, Fij , for a set of equiv- 
alent parameters q as follows [20( 



The Fisher matrix of q = (R, l a , ^h, ns) is the inverse of 
the covariance matrix of q (given in Tables I and II) . Note 
that the CMB shift parameters R and l a encode all the 
information on dark energy parameters. For any given 
dark energy model parameterized by the parameter set 
p^, the relevant Fisher matrix for p = (p x , ^de, 
uj m , ujb, ns) can be found using Eq. (|Bl[) . For the case 
most discussed in the literature, wx{z) = wq + w a (l — a), 
Px = ( w o, w a ). 

In order to find the Fisher matrix for (wq, w a , £Ide, 



2 We used a T C mb = 2.726, and R b = 30000C 6 /i 2 as defined in 
CAMB, noting that precision can be improved by updating these 
definitions. 
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f2fc, uj m , LOb, ns), the following derivatives are needed: 



OR = dT(zcMB) 
dwi dwi 
dhxR 1 




^k\ 1/2 r{z CM B) 



cosn[|O fc | 1 /2r(z CMB )] dTjzcMB) 

dQ.DE 

n k \ 1/2 

dT(z CMB ) 
OR 



wm[\n k \V*T(z C MB)] 
d\nR _ 1 1 
dftk 2f2 m 2£lfe 

COSnpflfcl^lXzcMB)] 

S inn[|fi fc |i/2r(z CMB )] 

9^6 



r(2CMs) 



0, 



9w m 
d In / a 



0. 



Oris 



d In Z Q 



= , fi ,1/2 COSn [Nfc| 1/2 r(^CMi3)] 

1 fc| sinn[|fi fe |i/2r(z CMB )] 
6>ln[ij" r 3 (2:c , MB)] 

cosnJO^T(£cMg)] 
sinn [|ilfc| 1 / 2 r(zcAfs)] 
51n[7J f s (^CMs)] 



= \n k \ 1/2 



dO.DE 

dT(zcMB) _ 
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dn k 



DE 



, n ,i/2 cosn[|n fc | 1 / 2 r(z C AfB)] 
1 fe| sinn [|n fc |Var(z CA , fl )i 



T(zcmb) 



20i 



91n[g r s (z CA f B )] 



<91nZ a 51n[iJ r s (2;(7Ms)] 



<9w m 
^ In l a 
duj b 
dl a 



dn s 



= 



dujb _ 



doj m 

d\n[H r s (z CMB )} 
duj b 



^ = (pi ^ w&) 

dpi 



9n 



= 1, 



s 



9ns 
dp; 



(pi ^ n s ), 



(B2) 



where Wi = (wo,w a ), and cosn(x) = cos(a;), x, cosh(x) 
for < 0, Ok = 0, and fife > respectively. 
Note that in the limit of O k — 0, 



dlni? = 91nr(z C MB) 
dO k 
din l a 

dO k dO k 
[T(zcmb)} 2 
6 



1 



do k 

dlnT(zcMB) 



inzcM B )f 

gln[g r s (z CM a)] 
dO k 



(B3) 



For the fiducial model considered in this paper, the 
Planck Fisher matrix for (wq, w a , Ode, O k , w m , w D , ns) 
is derived from the Planck covariance matrix of (i?, l a , 
Lu , ns) given in Table IV. 



APPENDIX C: USING OUR LIKELIHOOD FOR 
HAMILTONIAN MONTE CARLO 



While most cosmological codes use the standard 
Metropolis MCMC algorithm, there are other MC ap- 
proaches which may provide faster exploration especially 
in high dimensions. One example is Hamiltonian Monte 
Carlo (HMC) [U [H[ where each parameter 0i acquires 
a partner corresponding to a momentum variable iti , and 
the log-likelihood is regarded as a potential. The mo- 
menta are drawn from a univariate normal probability 
distribution and the next step in the MCMC exploration 
is chosen based on a Hamiltonian motion in this system, 
with total energy E = p 2 /2+x 2 {0) /2. At the end the mo- 
menta are marginalised over, which provides an ensemble 
of samples of the remaining parameters which is drawn 
from the posterior distribution. The main advantage of 
the HMC method is that the Hamiltonian motion nat- 
urally follows even complicated shapes of the posterior 
and in principle every proposal is accepted. The main 
drawback is that, in order to follow the trajectory, one 
needs to evaluate the gradient of the log-likelihood with 
respect to the parameters for dozens of steps, for every 
single proposal. Each proposal therefore requires hun- 
dred^) of likelihood evaluations if the gradient cannot 
be computed analytically. 

Since we have a reasonable approximation of the likeli- 
hood, we can instead use this approximation to compute 
the gradients. This means that the motion follows the 
(R,l a ,ujb,n B ) likelihood and at the end the approximate 
and the true likelihood are compared. If the true like- 
lihood is worse than the approximate one, then we can 
either assign the ratio as a weight to the new point (im- 
portance sampling) or test for rejection with the usual 
criterion (rejection sampling). If the true likelihood is 
better, then have to assign the ratio as a weight > 1. For 
this to work we must ensure that the approximate like- 
lihood does not exclude parameter regions that the true 
likelihood would allow. 

In our case we find that the procedure works quite well 
for the case where the Planck data is combined with the 
SN la and BAO data, since there the information loss 
is negligible 3 . Indeed, we find about 20% efficiency (ie 
roughly every 5th proposal is accepted, or correspond- 
ingly, the average weight of each point is 0.2), which 
is very good, especially since we can move a long dis- 
tance and obtain completely uncorrelated samples. Us- 
ing only the Planck data, we lose a lot of information, 
and less than 2% of the proposals are accepted. This is 
still not too bad, considering the complexity of the shape 
of the posterior, and that the resulting samples are com- 
pletely decorrelated. Additionally, burn-in is very quick 



3 We add additionally r and In A s to the set of our parameters, and 
augment the (R, l a , £^6, n B ) likelihood with their 2x2 covariance 
matrix. 
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for HMC and there is no need for initial runs to determine the optimal proposal matrix. 
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TABLE I: R, l a , Q b h 2 and n s estimated from Planck simulated data. 



Parameter mean rms variance 



R 



Q k 



1.7016 0.0055 

302.108 0.098 

0.02199 0.00017 

0.9602 0.0038 



TABLE II: Covariancc matrix for (R,l a ,Q b h 2 ,n s )h om Planck. 





R 


L 


Q b h 2 


n s 


Q k t^o 


R 


0.303492E-04 


0.297688E-03 


-0.545532E-06 


-0.175976E-04 


la 


0.297688E-03 


0.951881E-02 


-0.759752E-05 


-0.183814E-03 


Qbh 2 


-0.545532E-06 


-0.759752E-05 


0.279464E-07 


0.238882E-06 


n s 


-0.175976E-04 


-0.183814E-03 


0.238882E-06 


0.147219E-04 



TABLE III: Normalized covariance matrix for (7?, l a , Qbh 2 , n s )from Planck. 





R 


la 


Qbh 2 


n s 


Q k 


R 


1. 


0.553856 


-0.592359 


-0.832527 


la 


0.553856 


1. 


-0.465820 


-0.491026 


Q b h 2 


-0.592359 


-0.465820 


1. 


0.372425 


n s 


-0.832527 


-0.491026 


0.372425 


1. 



TABLE IV: Fisher matrix for (wo, w a , Qde, Qk, fc, UJb, ns) derived from the covariance matrix for (R, l a , Q b h 2 , n s ) from 
Planck. 





Wo 


W a 


n DE 


n k 




UJb 


ns 


wo 


.172276E+06 


.490320E+05 


.674392E+06 


-.208974E+07 


.325219E+07 


-.790504E+07 


-.549427E+05 


W a 


.490320E+05 


.139551E+05 


.191940E+06 


-.594767E+06 


.925615E+06 


-.224987E+07 


-.156374E+05 


^DE 


.674392E+06 


.191940E+06 


.263997E+07 


-.818048E+07 


.127310E+08 


-.309450E+08 


-.215078E+06 


Cl k 


-.208974E+07 


-.594767E+06 


-.818048E+07 


.253489E+08 


-.394501E+08 


.958892E+08 


.666335E+06 


UJm 


.325219E+07 


.925615E+06 


.127310E+08 


-.394501E+08 


.633564E+08 


-.147973E+09 


-.501247E+06 


UJ b 


-.790504E+07 


-.224987E+07 


-.309450E+08 


.958892E+08 


-.147973E+09 


.405079E+09 


.219009E+07 


ns 


-.549427E+05 


-.156374E+05 


-.215078E+06 


.666335E+06 


-.501247E+06 


.219009E+07 


.242767E+06 



